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We present the design and performance of a non-imaging concentrator for use in broad-band polarimetry 
at millimeter through submillimeter wavelengths. A rectangular geometry preserves the input polarization 
state as the concentrator couples // 2 incident optics to a 2tt sr detector. Measurements of the co-polar and 
cross-polar beams in both the few-mode and highly over-moded limits agree with a simple model based on 
mode truncation. The measured co-polar beam pattern is nearly independent of frequency in both linear 
polarizations. The cross-polar beam pattern is dominated by a uniform term corresponding to polarization 
efficiency 94%. After correcting for efficiency, the remaining cross-polar response is -18 dB. 
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1. Introduction 

The advent of bolometric detectors with background- 
limited sensitivity has important implications for astro- 
nomical instrumentation at millimeter and submillime- 
ter wavelengths. Once the detector phonon noise falls 
below the background from photon statistical fluctua- 
tions, further sensitivity gains can only be realized by 
collecting additional photons. A common implementa- 
tion increases the effective detecting area through an 
array of individual sensors, each coupled to the sky 
through an optical structure (feed horn, lenslet, phased 
antenna array) that restricts the system response to 
a single electromagnetic mode at the sensor. Single- 
moded systems achieve diffraction-limited angular res- 
olution with well-defined (Gaussian) beam profiles, but 
the large number of sensors required drives system-level 
complexity and cost. 

An alternative design uses a much smaller number of 
detectors capable of sensing multiple modes of the in- 
cident field. A multi-moded system uses non-imaging 
optics (a “light bucket” or similar feed structure) to fix 
the beam size independent of wavelength. Since the La- 
grange invariant is conserved, in order for the beam size 
to remain constant the number of modes must scale with 
wavelength, N mo ^ e = AQ/ A 2 , where A is the detector 
area, Ul is the solid angle, and A is the observing wave- 
length. The different electromagnetic modes form an 
orthogonal basis set and add incoherently at the detec- 
tor. The number of detectors required to reach a given 
sensitivity thus scales with etendu as (AM) -1 , inversely 
proportional to the number of modes. 

Multi-moded systems can be particularly useful for 
observations of the linear polarization of the cosmic mi- 


crowave background (CMB), where the signals of inter- 
est he well below the photon background noise. Po- 
larimetry adds an additional level of complexity to the 
optical system. Multi-moded systems typically require 
a concentrator or similar feed structure to couple light 
from the sky to the detector while conserving etendu. 
Polarization sensitivity can be achieved using a detec- 
tor sensitive to both polarizations, with a polarization 
diplexer {e.g. a wire grid analyzer) upstream of the de- 
tector. Since the polarization separation occurs before 
light enters the concentrator, the polarization properties 
of the concentrator do not matter. Alternatively, the 
detector itself may be sensitive to a single linear polar- 
ization. This avoids the need for a separate polarization 
diplexer, but requires the concentrator to preserve the 
incident polarization state without cross-polar mixing 
(see, e.g., the discussion in [1]). Instruments in this lat- 
ter class include the planned PIXIE [2], MUSE [3] and 
LSPE [4] missions. 

The design of non-imaging concentrators has been 
well studied for unpolarized observations [5-10]. Three- 
dimensional concentrators formed as a surface of revolu- 
tion can approach the theoretical limit on concentration, 
but the resulting azimuthal symmetry induces an unde- 
sirable cross-polar response. In this paper, we describe 
the design and performance of a non-imaging concen- 
trator using a rectangular geometry suitable for CMB 
polarimetry. 

2. Multi-moded rectangular concentrator 

Figure 1 illustrates the concentrator, derived from the 
requirements of the PIXIE instrument. We define a 
right-handed coordinate system [x,y,z\ with the de- 
tector centered in the xy plane and the concentrator 
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Fig. 1. Schematic of the off-axis non-imaging rectangular concentrator. The square entrance and exit apertures preserve 
incident polarization. Each of the 4 walls is an elliptical section to transform the // 2 incident beam to the 2 tt sr solid angle at 
the detector while preserving etendu. 


extending in the +z direction. Light from the // 2 
beam-forming optics is incident on the concentrator at 
a 15° angle from the z axis. Within these constraints, 
we use a rectangular geometry for the concentrator to 
minimize cross-polar mixing of the incident light. Each 


of the 4 walls of the concentrator is an elliptical sec- 
tion, optimized from an iterative ray-tracing algorithm 
to maximize the total throughput in the geometric optics 
limit. We define the ellipse in each plane as a quadratic 
equation 


Mi(x - ci) 2 + M 2 {z - c 2 ) 2 + 2 M 3 (x - q)(z - c 2 ) = 1 (top face) 

Mi(x — ci) 2 + M 2 (z — c 2 ) 2 + 2M%{x — ci) (z — c 2 ) = 1 (bottom face) 

Mi (y - ci) 2 + M 2 (z - c 2 ) 2 + 2 Ms(y - c 1 )(z - c 2 ) = 1 (side faces) (1) 


for the top, bottom, and side faces, respectively. Ta- 
ble 1 lists the coefficient values M and c. The two side 
panels (+y and —y) are identical, while the top and bot- 
tom panels (+£ and —x) have different elliptical sections 
to accommodate the off-axis design. The projected en- 
trance aperture thus appears approximately square. 

The exit aperture of the concentrator illuminates a 
pair of detectors, each sensitive to a single linear polar- 
ization. Each detector covers the full aperture; they are 
mounted one behind the other so that both are fully il- 
luminated. The two detectors are rotated by 90° such 
that one detector is sensitive to polarization in the x di- 
rection and the other is sensitive to polarization in the y 
direction. Each detector has etendu AQ = 4 cm 2 sr. We 
model the co-polar and cross-polar response of the con- 
centrator using a time-reversed Monte Carlo approach. 
For each linear polarization, we generate a set of 63 mil- 
lion outgoing rays, each originating from a random posi- 
tion on the detector. The ray angular distribution is uni- 
form in azimuth, but takes into account the additional 
cos 2 (6) angular response of a thin resistive detector at 


the concentrator exit aperture [1]. The distribution in 
polarization is uniform. A ray-tracing algorithm follows 
each ray through multiple bounces within the concen- 
trator (assumed to be perfectly reflective) to determine 
the orientation and polarization angle after leaving the 
concentrator. 

The PIXIE concentrator does not view the sky di- 
rectly, instead coupling to the beam-forming optics 
through an f/2 coupling mirror. We extend the ray- 
tracing algorithm to follow each ray out of the concen- 
trator and and through one reflection from the coupling 
mirror. The resulting set of rays approximates the beam 
pattern of the concentrator /coupling mirror combina- 
tion, which can then be used as input through the re- 
maining optical system. The co-polar and cross-polar 
response of any downstream beam-forming optics can 
readily be calculated, and and not treated here. 

Observations of the cosmic microwave background can 
span more than a decade in frequency, from 30 GHz to 
600 GHz. The corresponding number of modes within 
the concentrator ranges from N mo( i e = 4 at the low- 



3 


Table 1. Concentrator Elliptical Figure 


Face 

Cl 

C2 Ml 

m 2 

m 3 

Top 

-2.6043E+01 4.9987E+02 1.3455E-03 

1.0025E-04 

1.7423E-05 

Bottom 

3.5040E+02 

6.0139E+02 5.3181E-04 1.9821E-04 

-2.6392E-04 

Left 

-1.7033E+02 5.6460E+02 8.7479E-04 

1.3043E-04 

-2.1488E-04 

Right 

1.7033E+02 

5.6460E+02 8.7479E-04 

1.3043E-04 

2.1488E-04 


Note: Units are in mm. Values describe the larger scaled concentrator (§3) 


est frequency to 7V mo de = 1600 at the highest frequency. 
The concentrator thus samples the transition from a few- 
mode system to a highly over-moded system. Several au- 
thors have discussed methods to estimate the beam pat- 
tern for non-imaging concentrators in the few-mode limit 
[11-13]. Typically these methods require explicit mode 
matching. Although mode matching successfully recov- 
ers the full phase and amplitude information needed to 
generate far-field beam patterns, this method is com- 
putationally expensive and is difficult to apply for off- 
axis geometries where a full 3-dimensional treatment is 
required. To compare the modeled beam patterns to 
measurement, we employ a simpler approximation. We 
use the beam pattern from the ray-tracing algorithm to 
model the beam in the geometric optics limit. We then 
mathematically expand the beam using an orthogonal 
set of basis functions, and recompute the beam using a 
truncated subset of these functions to account for mode 
loss at low frequencies. The rectangular symmetry of 
the concentrator is well matched to the orthogonal basis 
functions of a two-dimensional Fourier transform. We 
thus bin the co-polar or cross-polar beam derived from 
the geometric optics model on a 6' x 6' grid, and Fourier 
transform the binned beam map to derive a set of com- 
plex Fourier coefficients. We sort the coefficients by an- 
gular frequency and retain only the lowest 7V moc ie values, 
setting the remaining modes to zero. We then Fourier 
transform back to real space to generate a model with 
the desired number of modes. Finally, we convolve the 
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Fig. 2. Block diagram for beam pattern measurement. 


mode-truncated pattern with the Airy pattern for the 
appropriate observing frequency to approximate the ef- 
fects of diffraction at the concentrator walls. 

3. Beam pattern measurements 

We measured the co-polar and cross-polar beam pat- 
terns of the rectangular concentrator in the Goddard 
Electromagnetic Anechoic Chamber (GEMAC). Figure 
2 shows the experimental setup. A transmitter within 
the GEMAC uses a standard rectangular gain feed to 
launch a single linear polarization toward a shaped re- 
flector. The reflector transforms the outgoing signal to a 
plane wave, effectively placing the transmitter at infin- 
ity. We include a coupling mirror in the GEMAC setup 
to convert the transmitted plane wave to the // 2 beam 
generated by the PIXIE coupling mirror, and measure 
the copolar and cross-polar beams from the concentra- 
tor/mirror system. The GEMAC coupling mirror differs 
from the PIXIE mirror in that the PIXIE mirror uses 
an elliptical surface to focus on the next mirror in the 
instrument while the GEMAC coupling mirror uses a 
parabolic surface to focus on the transmitter at infin- 
ity. Since the GEMAC setup excludes the remainder of 
the beam-forming optics, the measured beam patterns 
represent only the concentrator /coupling mirror combi- 
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Fig. 3. Schematic showing the concentrator and cou- 
pling mirror mounted for measurement in the GEMAC az- 
imuth/elevation table, showing the relation between the table 
and concentrator coordinate systems. 
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Fig. 4. Co-polar and cross-polar beam maps of the scaled concentrator / coupling mirror system. The measurement of the 
scaled concentrator at 90 GHz corresponds to 270 GHz for PIXIE, well within the geometric optics limit. The transmitter 
is polarized in the x direction for the copolar map and the y direction for the cross-polar map. The co-polar beam is well 
described by a tophat with diameter 14°. The cross-polar beam is dominated by a similar tophat at -12 dB (corresponding to 
polarization efficiency 0.94), with smaller-scale structure at the -18 dB level. 


nation and are not a measure of the final instrument 
beams on the sky. 

The concentrator /mirror system is mounted on a ta- 
ble and can be moved in azimuth and elevation. Figure 3 
show the relative orientation of the table and concentra- 
tor coordinate systems. Viewed from the detector in the 
concentrator xy plane, moving the table towards positive 
azimuth moves a source at infinity toward positive x (as 
viewed through the coupling mirror). Similarly, moving 
the table towards positive elevation moves a source at 
infinity toward positive y. To suppress reflected signals, 
all surfaces except the coupling mirror and concentrator 
aperture were covered with a single layer of Eccosorb 
AN72 microwave absorber. The absorber at the edges 
of the coupling mirror acts as a beam stop, truncating 
the corners to produce the desired tophat geometry. 

We modulate the transmitted power by electronically 
switching the transmitted power with a square wave 
at 20 Hz, and use a Thomas Keating THz absolute 
power meter [14] at the exit aperture of the concentra- 
tor to detect and synchronously demodulate the signal. 
The power meter is insensitive to polarization; a free- 
standing wire grid polarizer mounted between the power 
meter and the concentrator exit aperture provides polar- 
ization sensitivity. The polarizer is flush to the concen- 
trator exit aperture while the power meter is parallel to 
the exit aperture but 3 mm further back within a reflec- 
tive integrating cavity. The wire grid polarizer consists 
of copper-clad tungsten wires 40 microns in diameter 
(36 fi m tungsten and 1.3-2. 5 /am copper) on a 118 fi m 
pitch. With the grid wires oriented parallel to the ele- 
vation direction, the y polarization is reflected and the 
power meter is sensitive to x polarization. The grid can 
be rotated by 90° so that the wires are parallel to the 


azimuth direction. In this orientation, the x polarization 
is reflected and the power meter is sensitive to y polar- 
ization. The measured polarization isolation is better 
than -27 dB at frequencies 10 GHz to 300 GHz. 

With the transmitter fixed to broadcast x polariza- 
tion, we orient the polarizer grid parallel to the elevation 
direction to measure the XX co-polar response. We fix 
the elevation angle at —10° and continuously sweep in 
azimuth from —12° to +12°, recording data from —10° 
to +10° to avoid possible startup transients in the az- 
imuth drive. We then increment the elevation by 1° and 
repeat to map the concentrator /mirror response over a 
20° x 20° field. We then rotate the transmitter to the 
y polarization while leaving the polarizer grid parallel 
to the elevation direction to measure the YX cross-polar 
response. Finally, we rotate the wire grid polarizer by 
90° and repeat to measure the YY co-polar and XY 
cross-polar response. Since the off-axis concentrator is 
slightly different in the x and y directions, we expect 
slightly different responses for the 4 beam patterns. 

We wish to measure the beam patterns at frequen- 
cies representative of both the few-mode and highly- 
overmoded limits. The limited broadcast power avail- 


Table 2. Measured Map Power 


Frequency Max Co-Pol Max X-Pol Noise Co-Pol Noise X-Pol 


(GHz) 

(pW) 

(mW) 

(mW) 

(mW) 

11 

1560 

91 

1.4 

0.9 

29 

700 

73 

0.7 

0.5 

90 

210 

24 

0.4 

0.5 
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able at frequencies above 100 GHz makes direct beam 
mapping at these frequencies impractical. We thus in- 
crease the dimensions of the nominal PIXIE concentra- 
tor by a factor of three, and measure the larger (scaled) 
concentrator at lower frequencies where higher broad- 
cast power can be achieved. We measured all four beam 
patterns (XX, YX, YY, and XY) at each of frequencies 
11 GHz, 29 GHz, and 90 GHz, corresponding to CMB 
frequencies 33, 87, and 270 GHz. Table 2 summarizes 
the measurements. The maximum co-polar power mea- 
sured at each frequency ranged from 1500 /x W at 11 GHz 
to 210 y W at 90 GHz with rms noise 1.4 yW to 0.4 yW 
per 0.1° pixel. The noise floor is thus at approximately 
-30 dB at each frequency. 

At each frequency, we define the co-polar and cross- 


polar beam patterns as 



Cl 

o 

yT 

II 

P c (x,y)n 

(2) 

I Pc{x,y) dn 

G x (x, 7 /) — 

Px(x, y)Q 

(3) 

f P c (x,y) dfi 


where P c (x,y) is the co-polar power measured at each 
point (x,t/), P x (x,y) is the measured cross-polar power, 
and 


(fP c (x,p)dQ) 2 
f Pc (x,y) dn 


(4) 


is the co-polar beam solid angle. With this choice of 
normalization the beam solid angle is simply related to 
the beam pattern, 


Q = 



(5) 


but with a peak amplitude greater than unity, G max > 1 

[15]- 

Figure 4 shows the co-polar (XX) and cross-polar 
(XY) beam pattern measured at 90 GHz in the GEM AC, 
corresponding to 270 GHz for PIXIE. The co-polar beam 
is described by a tophat with diameter 14°. To lowest or- 
der, cross-polar beam resembles the co-polar beam but 
with amplitude reduced by 12 dB. Additional higher- 
order structure in the cross-polar beam is also visible, re- 
flecting the fourfold symmetry of the concentrator walls. 


Table 3. Beam Width 

Frequency Co-Pol Solid FWHM Equivalent 
(GHz) Map Angle (deg) Tophat 
(sr) (deg) 


11 

XX 

0.0502 

12.5 

14.5 

11 

YY 

0.0515 

13.4 

14.7 

29 

XX 

0.0483 

12.6 

14.2 

29 

YY 

0.0477 

12.4 

14.1 

90 

XX 

0.0470 

12.7 

14.0 

90 

YY 

0.0468 

12.7 

14.0 


Figure 5 compares the measured beam maps to the 
model. Each panel shows a one-dimensional slice 
through the beams along the x (azimuth) axis at y = 0. 
We show the co-polar XX beam pattern (transmitter 
broadcasting in x and concentrator receiving in x) as 
well as the cross-polar YX pattern (transmitter broad- 
casting in y and concentrator receiving in x). The co- 
polar beams agree well with the model and show the 
expected broadening from mode loss at low frequencies. 
The cross-polar response is reduced by ~12 dB from the 
co-polar response and is dominated by the same tophat 
structure as the co-polar beam. We observe similar re- 
sults for slices in elevation (x = 0) (Figure 6). Note 
that the right-most panels in Figures 5 and 6 corre- 
spond to horizontal or vertical slices through the 90 GHz 
beam shown in Figure 4. The measured beam patterns 
show the qualtitative features predicted by the mode- 
truncated model. Differences between the measured and 
modeled beam patterns are dominated by an amplitude 
scale factor which increases at lower frequency. Reflec- 
tions from the mounting plate beneath the concentra- 
tor/mirror system (Figure 3) could account for such an 
effect. 

Diffraction and mode truncation combine to apodize 
the beam edges. The full width at half maximum 
(FWHM) of the measured beam is 12.7°. Power out- 
side the half-power point broadens the beam so that the 
FWHM does not fully describe the co-polar beam. Ta- 
ble 3 shows the FWHM for each co-polar beam, along 
with the diameter of an ideal circular tophat with the 
same solid angle as the measured beam. Although the 
FWHM is nearly independent of frequency, the equiva- 
lent tophat diameter increases slightly at the lowest fre- 
quency, showing the effect of mode loss and diffraction. 


3. A. Polarization Efficiency and Cross-Polar Residu- 
als 

A cross-polar response proportional to the co-polar 
beam is equivalent to a loss in polarization efficiency, 
but does not generate systematic error when mapping 
the polarized sky. Figure 7 illustrates this point. With 
the concentrator held stationary, we record the received 
power at 11 GHz while the transmitter rotates in po- 
larization angle 0. The measured power shows the ex- 
pected sin(2$) modulation with a polarization efficiency 
94% corresponding to an on-axis cross-polar response of 
-12 dB. 

We estimate the non-uniform cross-polar response by 
fitting each cross-polar beam map to the corresponding 
co-polar map, 


JG X G C dtt 
a ~ Jofdn 

and computing the rms amplitude of the residual 


(6) 


Rx(x, y) = G x (x, y) - aG c (x, y ) . (7) 


Figure 8 compares the residual cross-polar response to 
the co-polar beam measured at 90 GHz. The beams 
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Fig. 5. Measured beam patterns for the scaled concentrator /coupling mirror system compared to model predictions in the x 
direction. The solid black line shows the co-polar (XX) beam while the dotted line shows the cross-polar (YX) beam. Data are 
taken with y fixed at 0°. Grey lines show the model predictions after accounting for diffraction and mode loss. Measurements 
at 11, 29, and 90 GHz correspond to 5, 34, and 324 modes respectively, spanning the transition from a few-mode system to a 
highly over-moded system. 
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Fig. 6. Measured and model beam patterns for the copolar (XX) and cross-polar (YX) beams in the y direction. Data are 
taken with x fixed at 0°. Model and measured data are as in Fig 5. 
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Fig. 7. On-axis response vs polarization angle measured at 
11 GHz. 


are shown on a linear scale to highlight the cross-polar 
structure. The residual cross-polar response shows the 
effects of the square aperture and walls. The asymmetry 
in the x direction results from the corresponding asym- 
metry in the concentrator due to the off-axis design (Fig- 
ure 1). Table 4 lists the polarization efficiency and rms 
cross-polar residuals measured for both polarizations at 
all three frequencies. Averaged over the full beam, the 
polarization efficiency is 94% and the rms cross-polar 
residual is -18 dB. 

4. Conclusion 

We have designed and tested an off-axis rectangu- 
lar non-imaging concentrator for use with multi-moded 
polarization-sensitive optical systems. The concentrator 
preserves linear polarization and conserves the optical 
etendu while transforming the beam from // 2 at the 
entrance aperture to 2tt sr at the detector (// 0.5). The 
measured co-polar beam pattern is nearly independent 
of frequency in both linear polarizations. The co-polar 
beam is well described by a tophat of angular diameter 
14° with minimal internal structure on smaller angular 
scales. The cross-polar response is dominated by a uni- 
form tophat, leading to polarization efficiency 94%. Af- 


Table 4. Polarization Properties 


Frequency Co-Pol Cross-Pol Polarization Cross-Pol 
(GHz) Map Map Efficiency (dB) 


11 

XX 

YX 

0.94 

-19 

11 

YY 

XY 

0.95 

-20 

29 

XX 

YX 

0.93 

-17 

29 

YY 

XY 

0.93 

-15 

90 

XX 

YX 

0.97 

-20 

90 

YY 

XY 

0.94 

-17 
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Fig. 8. Co-polar beam and cross-polar residual measured at 
90 GHz. Contours show the XX co-polar beam on a linear 
scale at amplitudes 0.01, 0.05, 0.1, 0.5, and 1. The YX cross- 
polar residuals are shown in color, with red positive and blue 
negative. 


ter removing the uniform term, the residual cross-polar 
response is typically -18 dB. 

Several improvements are possible. The design pro- 
cess restricted the figure for each wall of the concen- 
trator to remain an ellipse, while varying the three foci 
to optimize the total throughput (co-polar and cross- 
polar). A more ambitious design could add additional 
free parameters to the figure for each wall, changing the 
simple ellipse to (e.g.) a low-order polynomial. Simi- 
larly, the concentrator could be optimized to maximize 
the co-polar throughput while simultaneously minimiz- 
ing the cross-polar response. We have not yet attempted 
such improvements. 

We measured the co-polar and cross-polar beam pat- 
terns for an optical system consisting of the concentra- 
tor plus a coupling mirror. The results are thus rep- 
resentative of the beam patterns expected at output of 
the PIXIE instrument and do not represent the PIXIE 
beams on the sky. The full PIXIE optical system in- 
cludes additional beam-forming mirrors to reduce the 
tophat diameter to the design value 2.6° [2]. The beam- 
forming optics are optimized to further symmetrize the 
beam patterns. 
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